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Abstract 

The / and d electron density of states of the one-dimensional Falicov- 
Kimball model are studied in the weak-coupling limit by exact diagonalization 
calculations. The resultant behaviors are used to examine the d-electron gap 
(A(i), the /-electron gap (Aj), and the /d-electron gap (Aj^) as functions 
of the /-level energy Ef and hybridization V. It is shown that the spinless 
Falicov-Kimball model behaves fully differently for zero and finite hybridiza- 
tion between / and d states. At zero hybridization the energy gaps do not 
coincide (A^^ / Aj ^ ^fd)i and the activation gap Aj^ vanishes discontin- 
uously at some critical value of the /-level energy Efc- On the other hand, 
at finite hybridization all energy gaps coincide and vanish continuously at the 
insulator-metal transition point Ef = Efc- The importance of these results for 
a description of real materials is discussed. 

PACS nrs.:75.10.Lp, 71.27.+a, 71.28.+d, 71.30.+h 
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1 Introduction 



The Falicov-Kimball model (FKM) has become, since its introduction in 1969, one 
of the most popular examples of a system of interacting electrons with short-range 
interactions. The model was originally proposed to describe metal-insulator transi- 
tions [1] and has since been investigated in connection with a variety of problems 
such as binary alloys [2], the formation of ionic crystals [3], and ordering in mixed- 
valence systems [4] . It is the latter language we shall use here, considering a system 
of localized / electrons and itinerant d electrons coupled via the local interaction U 
and hybridization V . The Hamiltonian of the spinless FKM is 

H = Y: Ujdfdj + UY: Ufidtdi + EfY: Ufi + ^ E dtfi + h.c, (1) 

ij i i i 

where fj^, fi arc the creation and annihilation operators for an electron in the local- 
ized state at lattice site i with binding energy Ef and d'l, d,i arc the creation and 
annihilation operators of the itinerant spinless electrons in the ci-band Wannier state 
at site i. 

The first term of (1) is the kinetic energy corresponding to quantum-mechanical 
hopping of the itinerant d electrons between sites i and j. These intersite hopping 
transitions are described by the matrix elements tij, which are —t if i and j are the 
nearest neighbors and zero otherwise (in the following all parameters are measured 
in units of t). The second term represents the on-site Coulomb interaction between 
the (i-band electrons with density Ud = \ dfdi and the localized / electrons with 
density Uf = \ Y.iftfii where L is the number of lattice sites. The third term stands 
for the localized / electrons whose sharp energy level is Ef. The last term represents 
the hybridization between the itinerant and localized states. 

In this paper we perform exhaustive numerical studies of the spectral properties 
of the FKM with and without hybridization. While the static properties of the FKM 
are well understood at present [5, 6] (including the picture of valence transitions), 
the dynamical properties of the model are still unclear. Even, the spectral properties 
of the / electrons are not understood satisfactorily nor for y = 0, where only a few 
exact results are known for the infinite-dimensional systems [7, 8]. No exact results 
are known for nonzero hybridization and T = 0, with the exception of numerical 
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results obtained on very small clusters [9]. The first systematic study of dynamical 
properties of the FKM with hybridization has been performed recently by Craco [10] 
within so called static approximation. He studied the model in the strong-coupling 
limit (large U) and found that the system is insulating for small values of hybridiza- 
tion (the size of the / and rf-electron gap coincides, including the case V = 0) and 
with increasing V an insulator-metal transition takes place at some critical value 
of V — Vc- In our preceding paper we reexamined these strong-coupling results by 
exact-diagonalization calculations and obtained fully different conclusions [11]. In 
particular, we have found that for zero hybridization the gaps in the / and d electron 
density of states do not coincide, and almost all / electron spectral weight is located 
outside the d electron subbands. For nonzero hybridization the / and d electron 
gaps coincide, for both the symmetric and unsymmetric case, but no insulator-metal 
transition driven by hybridization is observed in one as well as in two dimensions. In 
this contribution we extend numerical calculations to the opposite limit (the weak- 
coupling limit). From this point of view the paper represents the first systematic 
(exact) study of dynamical properties of the FKM in the weak-coupling limit. Here, 
the special attention is devoted to the behavior of the / and ci-electron density of 
states with increasing Ef. Such an analysis is very important since a parametrization 
of Ef with applied pressure p [12] can, in principle, provide an interpretation of some 
experimental data, e.g., the behavior of the activation gap (the gap between the occu- 
pied and unoccupied states) with increasing p. In our previous paper [6] the problem 
of pressure dependence of the activation gap was analysed through the behavior of 
the energy gap in the rf-electron spectrum of the FKM and many similarities with 
the experimental data on the valence fluctuating compound SmBg were found. These 
similarities are however only qualitative, since in the correct analysis one should take 
into account also the behavior of the /-electron spectral functions. This was done 
by Park and Hong [9] and a nice correspondence of theoretical and experimental re- 
sults has been obtained. Unfortunately, these results have been obtained on a very 
small cluster (consisting of only eight sites) and thus cannot be considered as definite. 
Here we reexamine the behavior of the activation gap in the FKM for both V — 
and V > 0. At zero hybridization we are able to present results on relatively large 
clusters for both the d and / electron density of states, and thus our results can be 
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extrapolated satisfactorily to the thermodynamic limit. For nonzero hybridization 
we were able to perform small-cluster exact-diagonalization calculations on lattices 
only slightly larger than used by Park and Hong [9] (L ~ 12), however a fundamental 
different behavior of the model is observed already on such small clusters. 

2 Results and discussion 

Let us start the discussion of our results with the case V = 0. In Fig. 1 we present 
numerical results for the d and / electron density of states obtained for U = 0.6 
and several values of the /-level energy Ef. Since the ground state configurations of 
the FKM in the weak coupling limit are well-known [6, 13] (the most homogeneous 
(insulating) configurations for \Ef \ < Efc ~ 1.34, and the phase separated (metallic) 
configurations for \Ef\ > Efc), the d-electron density of states can be calculated 
directly from the single particle spectrum of the FKM model for y = 0. For the 
periodic configurations this can be done analytically in the thermodynamic limit [14], 
and for an arbitrary /-electron concentration a numerical diagonalization is possible 
on very large clusters (L = 64000). Unhke this case, the /-electron density of states 
has to be calculated by exact diagonalization calculations (Lanczos method [15]), 
thereby the cluster sizes are strongly limited (L < 24). To minimize the finite-size 
effects on the /-electron spectra we have performed numerical calculations for several 
cluster sizes at each value of Ef. The typical behaviors are shown in Fig. 2 for two 
selected values of Ef. It is seen that finite-size effects are small and thus already 
results obtained for L = 24 can be used satisfactorily to represent the behavior of 
large systems. Although we have used in the numerical calculations the relatively 
large value of e = 0.01 for the resolution [16], the formation of the gap in the /- 
electron density of states at the Fermi level (in all examined cases the Fermi level is 
located between the first and the second peak) is apparent. It is also apparent that 
the gaps in the d and / electron spectra do not coincide. Moreover, one can see (in 
accordance with results obtained in the strong-coupling limit [11]) that practically all 
/-electron spectral weight is located inside the principal (i-electron gap (the gap at the 
Fermi energy). With increasing Ef both the principal gap as well as the /-electron 
spectrum shift to higher energies. Since the finite size effects on the /-electron spectra 
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are negligible for the cluster sizes used in our numerical calculations (L = 24), and 
even the d-electron spectra can be obtained exactly on much larger clusters, one can 
try to construct the Ef dependence of all relevant energy gaps. In particular, we have 
calculated the d-electron gap (the gap in the d-electron spectrum at the Fermi 
energy), the /-electron gap Af (the gap in the /-electron spectrum at the Fermi 
energy), and the /rf-electron gap Afd (the gap between the occupied / subband and 
the empty d subband (the activation gap). To obtain the ci-electron gap it is 
sufficient to know the /-electron distribution that minimizes the ground-state energy 
of the FKM for given Ef and U. As mentioned above, the ground-state configurations 
of the FKM in the weak-coupling limit are the most homogeneous configurations for 
< Efc, and the phase separated configurations for \Ef \ > Efc- While the most 
homogeneous distribution of / electrons can be easily generated for arbitrary Ef{nf), 
the /-electron distribution in the phase separated configuration has to be determined 
numerically. This distribution can be found in principle exactly, as shown in [6] or 
approximately as done in [13]. Here we adopt the latter method, since it allows 
to treat several times larger clusters (L ~ 300) and still to keep the high accuracy 
of computations. The results of numerical calculations for A^ obtained using this 
method on the cluster consisting of 240 sites are presented in Fig. 3. It is seen that 
Ef dependence of the d electron gap qualitatively mimics the pressure dependence 
of the activation gap in SmBg [17] (we note a parametrization of Ef with pressure) . 
With increasing pressure (the /-level position Ef) the gap decreases and vanishes 
discontinuously at some critical pressure pc (Efc). The finite-size scaling analysis 
that we have performed for a wide range of L values {L = 40, 60 . . . 240) showed 
that Ad is practically independent of L for all \Ef \ < Ef^. On the other hand the 
same analysis performed for \Ef\ > Efc revealed that the finite-size effects are still 
present in this region, but a disappearance of A^ can be satisfactorily confirmed by 
extrapolation of results obtained for different L. 

To determine the Ef dependence of the /-electron gap Af , and the /d-electron 
gap Afd we have performed exhaustive numerical studies of the /-electron density 
of states for a wide range of Ef values. In particular, we have calculated the /- 
electron density of states for each Ef from the interval [0,1.4] with the step AEf — 
0.05, and with the step AEf — 0.01 for Ef near the point of the insulator-metal 
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transition. To reveal the finite-size effects on A/^ and A/ the calculations have 
been done for several cluster-sizes (L = 12, 16, 20, 24) at each Ef. It was found 
that finite-size effects on A/ and A/^ are negligible for \Ef\ < Efc, while small 
finite-size effects have been observed for Ef > Efc- The resultant Ef dependence of 
Af and Af^ is presented in Fig. 3. It is seen that Af^ exhibits the same behavior 
as Ad, of course with one exception and namely that A^/ ~ Ad/2. Since Afd is 
the gap between the occupied (/) and unoccupied (d) states (the activation gap) 
its behavior can be directly compared with the behavior of the activation gap in 
SmBe [17]. Although data for Afd are more scattered, both activation gaps exhibit 
quahtatively the same behavior. This indicates that the spinless FKM, probably the 
simplest model of correlated electrons can, in principle, provide the correct physics 
for describing properties of real materials. 

Let us now examine what happens if the hybridization is switched on. One could 
expect that just the hybridization will improve the accordance between the theoretical 
and experimental results for the activation gaps, since the hybridization generally 
smears behaviors. This conjecture support also results by Park and Hong [9] obtained 
for the FKM with hybridization on a very small cluster consisting of only eight 
sites. Their results obtained for the Ef dependence of the activation gap and the d- 
electron concentration Ud are summarized in Fig. 4. It is seen a nice correspondence 
of theoretical and experimental data for both the activation gap as well as the d- 
electron concentration. Unfortunately, these results have been obtained on a very 
small cluster and so cannot be considered as definite. Even, the actual behavior of 
the activation gap and the ri-clectron concentration Ud on large lattices can be fully 
different from one presented in Fig. 4. Indeed, a comparison of their results obtained 
for Ud d.t V = with our results presented in Fig. 3b reveals fully different behavior 
of Ud- While Ud calculated by Park and Hong [9] for L = 8 is constant for all values 
of Ef from Ef — 0.8 to the insulator-metal transition point Efc — 1.34, our results 
calculated on much larger clusters (L = 240) exhibit strong dependence on Ef. In 
the next we reexamine exactly the behavior of the FKM with hybridization on the 
cluster consisting of eight sites, as well as on clusters slightly larger (L = 10, 12). We 
present new results that significantly improve results obtained by Park and Hong [9] 

To show the hybridization effects on Ad,Af and A^^ we have calculated exactly 
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the / and d electron density of states for y = and V — 0.02 on a small finite cluster 
consisting of L = 8 sites. We chose this cluster size to be compatible with results 
of Park and Hong [9]. The resultant behaviors are presented in Fig. 5. The most 
prominent difference between the results obtained for F = and V = 0.02 is that 
a nonzero spectral weight appeared in the (/-electron density of states at the Fermi 
level for finite hybridization. Of course, this fact will change dramatically the picture 
discussed above for V = 0. This is clearly demonstrates in Fig. 6a, where the Ef 
dependence of A^, A/ and A fa is displayed for V — 0.02 and L — 8. It is seen that 
all gaps coincide for nonzero V, what strongly contradicts to the case oiV — 0. For 
a comparison we have displayed in Fig. 6a also the behavior of the single particle 
excitation energy defined as A^ = Eg{L + 1) + Eg{L — 1) — 2Eg{L), where Eg{N) 
is the ground state energy for A^" electrons. As one could expect (on the base of the 
strong-coupling results [11]) the single particle excitation energy A^ coincides with 
Ad, Af and Afd for nonzero hybridization. This result is very important from the 
numerical point of view since the single particle excitation energy can be calculated 
easily by other methods (that allow to treat much larger clusters, e.g., the density 
matrix renormalization group method (DMRG)), and so it can be used satisfactorily 
for describing conducting properties of the model. As one can see from Fig. 6a the 
conducting properties of the FKM with hybridization are described very inaccurate 
within the small-cluster exact-diagonalization calculations. In the region where the 
metallic phase has been identified for V — 0, the single particle excitation energy in- 
creases and it is not clear if this is a consequence of finite hybridization or finite-size 
effects. Since calculations performed for L = 10 and L = 12 revealed relatively large 
finite-size effects we have decided to use DMRG method to verify the actual behavior 
of Ag. In comparison to exact diagonalization calculations the DMRG method allows 
to treat several times larger clusters (L ~ 50) and still to keep the high accuracy of 
computations. We typically keep up to 128 states per block, although in the numer- 
ically more difficult cases, where the DMRG results converge slower, we keep up to 
400 states. Truncation errors [18], given by the sum of the density matrix eigenvalues 
of the discarded states, vary from 10~^ in the worse cases to zero in the best cases. 
The typical behavior of A^ is displayed in Fig. 6b for two finite clusters of L = 32 
and L = 36 sites. Although the finite-size effects are still present for cluster sizes 
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treated by DMRG method the metaUic region seems to be satisfactorily verified also 
for nonzero hybridization. However, there is one important difference between the 
case V — and the case V — 0.02. While for y = both A,^ and Afd vanish discon- 
tinuously, they seem to vanish continuously for nonzero hybridization. This indicates 
that the spinless FKM behaves fully differently for zero and non-zero hybridization 
and fully different are also the corresponding pictures of insulator-metal transitions. 
Since rare-earth compounds exhibit both types of insulator-metal transitions (dis- 
continuous as well as continuous) their different behavior under the external pressure 
can be interpreted directly as a consequence of absence (presence) of hybridization 
in a given material. According to this conjecture, e.g., SmBe that exhibits a discon- 
tinuous pressure-induced insulator-metal transition should be a material with zero 
hybridization between d and / states. Indeed, SmBe is the highest symmetry of the 
Oh point group and the on-site hybridization between d and / states is forbidden by 
inversion symmetry [19]. This confirms that the spinless FKM, in spite of its rela- 
tively simplicity, can yield the correct physics for a description of real materials, e.g. 
SmBg. Besides of a qualitative correspondence between the theoretical and experi- 
mental data for the activation gap, we can state another example that supports this 
conclusion. In Fig. 7 we have displayed in detail the d and / electron density of states 
near the Fermi level for several values of Ef {Ef — 0, 0.1, 0.2). For this set of Ef val- 
ues the ground state is characterized hy Uf — 0.5 what models very well the situation 
in SmBg at low temperatures [6] . This electronic structure consisting of two wide and 
two narrow subbands strongly mimics the electronic structure that we have proposed 
several years ago for a description of low temperature resistivity data of SmBg [20]. 
The analysis of these data showed that they can be explained by introducing a fine 
structure (consisting of two narrow bands) into the principal gap, in accordance with 
our results presented in Fig. 7. Of course, all these comparisons are only qualitative 
since our results have been obtained for the one-dimensional system at T = 0, while 
the real experimental systems are three dimensional and measurements are done at 
finite temperatures. In future work, we plan to perform a similar analysis in higher 
dimensions and T 7^ 0. Moreover, we also want to examine the infiuence of other 
factors (the electron-phonon interaction, the orbital dynamics, etc.) that have been 
neglected in this version of the model. 
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In summary, the / and d electron density of states of the spinless Fahcov- Kimball 
model were studied in the weak-coupling limit by exact diagonalization calculations. 
The resultant behaviors were used to examine the d-electron gap A,^, the /-electron 
gap Af, and the /^-electron gap Afd as functions of the /-level energy Ef and 
hybridization V. It was shown that the spinless Falicov-Kimball model behaves 
fully differently for zero and finite hybridization between / and d states. At zero 
hybridization the energy gaps do not coincide (A^ ^ Af ^ A/^), and the activation 
gap A/d vanishes discontinuously at some critical value of the /-level energy Ejc- 
In all examined cases practically all /-electron spectral weight is located outside 
the d-electron subbands. On the other hand, at finite hybridization all energy gaps 
coincide and vanish continuously at the insulator-metal transition point Ej — Ejc- 
The results obtained are compared with experimental behaviors observed in SmBg. 
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2/7021/20 and the Science and Technology Assistance Agency under Grant APVT- 
51-021602. Numerical results were obtained using the PC- Farm of the Slovak Academy 
of Sciences. 
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Figure Captions 

Fig. 1. The / and d electron density of states of the FKM obtained for V — 
0, [/ = 0.6 and several different values ol Ef. The results for the / electron density of 
states are plotted for L — 24, while the results for the d electron density of states are 
plotted for L — 64000. The corresponding /-electron concentrations are: n/ = 0.5 for 

= Q,nf = 4/10 for Ef = 0.5, n/ = 1/3 for Ef = 0.8, and n/ = 1/4 for Ef = 1.2. 

Fig. 2. The / electron density of states of the FKM obtained for F = 0, C/ = 0.6 
and two different values of L and Ef. 

Fig. 3. The d-electron gap A^, the /-electron gap A/, the /(i-electron gap A/^, 
and the (/-electron concentration Ud as functions of the /-level energy Ef calculated 
for V = and U = 0.6. The solid lines are guides to the eye. 

Fig. 4. The activation gap Afd and the d-electron concentration Ud as functions 
of the /-level energy Ef calculated for U = 0.6 and L = 8 (Ref. 7). 

Fig. 5. The / and d electron density of states of the FKM obtained for U = 
0.6, Ef = 1.2, L = 8 and two different values of V. The arrow denotes the most 
prominent difference between the results obtained for V = and V = 0.02. 

Fig. 6. a) The rf-electron gap A^^, the /-electron gap Aj, the /rf-electron gap Afd, 
and the single particle excitation energy Ag as functions of the /-level energy Ef 
calculated for V = 0.02, L = 8 and U = 0.6. b) The single particle excitation energy 
Ag as a function of the /-level energy Ef calculated by DMRG method for the same 
values of V and U. The sohd and dashed lines are guides to the eye. 

Fig. 7. The / and d electron density of states of the FKM obtained for V — 
0,U — 0.6 and several different values of Ef. The results for the / electron density 
of states are plotted for L — 24, while ones for the d electron density of states are 
plotted for L — 64000. In all cases the /-electron concentration is equal to Uf — 0.5. 
To see clearly the /-electron gap, the /-electron density of states has been plotted 
with high resolution (e = 0.001). 
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